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^.^ Abstract 

rS In light of the recent discovery of an approximately 126 GeV Higgs 

I boson at the LHC, the particle physics community is beginning to explore 

fS . the possibilities for a next-generation Higgs factory particle accelerator. 

U In this report we study the s-channel resonant Higgs boson production 

i-C and Standard Model backgrounds at a proposed M^p" collider Higgs fac- 
tory operating at center-of-mass energy y^i = Mh with a beam width of 

,__! 4.2 MeV. We study PYTHIA-generated Standard Model Higgs and back- 

^ ground events at the generator level to identify and evaluate important 

f^ channels for discovery and measurement of the Higgs mass, width, and 

t~^ branching ratios. We find that the H° — > bb and H^' — >• WW* channels 

Cn are the most useful for locating the Higgs peak. With an integrated lumi- 

'•'i nosity of 1 fb~^ we can measure a 126 GeV Standard Model Higgs mass 

^h accurately to within 0.25 MeV and its total width to within 0.45 MeV. 

C^ Our results demonstrate the value of the high Higgs cross section and 

CO narrow beam resolution potentially achievable at a muon collider. 
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1 Introduction 

The Higgs field and its associated resonance, the Higgs boson were added to the 
Standard Model to solve the problem of how bosons and leptons acquire mass 
without breaking gauge invariance. The Higgs field and mechanism have since 
become relevant to fundamental theories about dark energy, dark matter and 
vacuum stability. The Higgs boson has been sought after for over forty years 
and until recently was the only remaining undiscovered fundamental particle 
predicted by the Standard Model. 

In July 2012, the CMS and ATLAS experiments at the Large Hadron Collider 
in CERN announced the discovery of a new particle consistent with Standard 
Model predictions of a Higgs boson with a mass between 125 and 126 GeV. While 
further measurements have continued to indicate that the particle is a Higgs 
boson, it remains to be seen whether all of its properties are consistent with the 
Standard Model. [7] The LHC will not be able to make the precise measurements 
needed to test the Standard Model to its limits. [3 The Standard Model predicts 
a width of F^ = 4.21+3g|MeF for a Higgs particle with a mass of 126 GeV 
but the LHC will be limited to a mass measurement with uncertainty on the 
order of tens of MeV due to detector resolution. [2] [3] Some decay channels, such 
as H'^ — > bb have very high irreducible physics backgrounds at the LHC, making 
branching fractions and partial widths very difficult to measure. [6j[3J 

The fundamental importance of the Higgs boson and its exact properties 
motivates the construction of a new particle accelerator 'Higgs factory'. There 
are many proposed candidates for a next-generation Higgs factory, including 
linear and circular e~^e~ colliders, 77 colliders, plasma wakefield accelerators 



and a circular /i^/^^ collider. [T] This report examines the potential ability of 
a proposed muon collider to fill that role and probe the Standard Model to 
its limits. We examine the physics backgrounds relevant to a muon collider 
operating at the Higgs s-channel resonance and explore an energy scanning 
search strategy for locating the narrow Higgs peak. We find that the high beam 
energy resolution and ability to use s-channel resonance Higgs production at a 
muon collider make it an attractive option for further research and development. 
In this report we assume a Standard Model Higgs. Unless stated otherwise, 
we assume a mass and width of 

Mh = 126.0Gey, Th = 4.21MeF (1) 

1.1 S-Channel Resonant Higgs Boson Production 

The Higgs boson's resonant production cross section is given by the Breit- 
Wigner formula. For a center of mass energy Vs, this is given by [S] 

{s-Mlf+T^Ml 

Vh is referred to as the 'width' of the Higgs peak. The peak value of the cross 
section, using Standard Model values for the width and branching fractions of 
a 126 GeV Higgs, is approximately QAph. The observable cross section is in 
practice the convolution of this Higgs peak with the energy distribution of the 
collider. We assume that the distribution of the center of mass energy is a 
Gaussian and unless otherwise stated, use a beam with a standard deviation in 
•\/I of 4.2 MeV, roughly the same as the Higgs peak. As will be shown later, 
this is an optimal width for discovering the Higgs. To calculate cross sections 
and to fit simulated data we numerically convolute the Higgs Breit-Wigner with 
a Gaussian. The peak value of the smeared cross section is 28.3p6. 

2 Muon Collider as a Higgs Factory 

2.1 Lepton Mass Coupling 

The Higgs mechanism couples to the square of a leptons's mass so s-channel 
resonance Higgs production is enhanced by a factor of 4.28 x 10'* in a muon 
collider as compared to an electron collider [4]. 

2 

9J1}^ oc ^ = 4.28 X 10^ m (3) 

dHee rn^ 

In e+e~ colliders the only feasible channels for Higgs production are Higgs- 
strahlung {e^e~ — ?> ZH) and vector boson fusion (e^e~ — ^ Hv^Ve), which have 
lower cross sections, higher physics backgrounds and do not allow for direct 
measurement of the Higgs mass and width. 

2.2 Beam Energy Resolution 

Many properties of a muon collider make it an attractive option for a Higgs 
factory. The high mass of the muon compared to the electron, a ratio of about 



200, suppresses radiative effects like beamstrahlung, allowing for a much greater 
beam energy resolution than an electron-positron collider. Estimates for achiev- 
able beam energy resolution at center of mass energies near 126 GeV are on the 
order of 1 to 10 MeV, surpassing the capabilities of any e^e^ machine by a 
factor of ten or greater. [^ This high energy resolution allows for direct mea- 
surement of the Higgs width by 'scanning' the beam energy across the s-channel 
resonance peak. 

With a beam energy resolution equal to the Higgs total width, the Higgs 
production cross section is a{^~^fi~ — >■ iJ'^) = 28.3p6.[4j Thus it is feasible to 
directly measure the Higgs width by 'scanning' the beam energy across the Higgs 
peak. Estimates for the luminosity of a 63 GeV per beam muon collider are on 
the order of L = 10^^cm~^s~^, suggesting the production of about 50,000 to 
100,000 Higgs bosons per year when taking data at the s-channel resonance. 
This would be competitive with a 126 GeV per beam linear e~^e~ collider with 
L = 2.0 X 10^^cm~^s~"'^ producing Higgs bosons in the Higgs-strahlung and 
vector boson fusion processes. [5] 

3 Machine Induced Background 

One disadvantage of a muon collider is that a large number of muons decay 
in the beamline, causing significant background. While we are mostly ignoring 
this for the purposes of this report, we do include a cos{9) restriction on our 
particles corresponding to a cone angle of 20 degrees. This is a conservative 
estimate motivated by the need for the cone to be thick enough to provide some 
shielding. Figure [l] is a conceptual rendering of MCDRCalOl, a detector with 
a dual-readout capability that measures both ionizing radiation and Cerenkov 
photons in the totally active calorimeters. This detector will be used in future 
simulation studies but is not used as a reference for this report. Simulated Higgs 



event displays are provided in Appendix A. 3 



4 Physics Background 

The most significant background for s-channel resonance Higgs production at a 
muon collider is the production of Z bosons. The Higgs cross section, smeared 
by a 4.2 MeV beam is 28.3 pb. The cross section of the Z background is 376 
pb, but 20.04% of these Z's decay into pairs of neutrinos and a photon, bringing 
the cross section to 301.4 pb and S/\/~B to 1.63. This cross section remains 
essentially flat in the region around the Higgs peak and will be treated as such 
in this report. Figure [2] shows simulated data of a scan across a 126.0 GeV 
Higgs peak counting all events except for Z^ — >■ vivg. The data is fitted to 
a Breit-Wigner convoluted with a Gaussian with three free parameters; Vh: 
Mh and Br{H^ — > X). The fixed parameters are the background cross section 
a{Z^ — >• X), the beam width ateam and the total integrated luminosity C. The 
fit gives a width of 4.56 ± 1.52 MeV, an error in the mass measurement of 
0.13 ± 0.16 MeV and a branching ratio of 0.96 ± 0.04. 




Figure 1: Conceptual rendering of MCDRCalOl detector. The silicon tracker 
(orange) outer radius is 122 cm. The electromagnetic and hadronic calorimeters 
(gray) are made of Bismuth Germanium Oxide, or BGO. The EM calorimeter 
has an inner radius of 125 cm and has ten 2cm thick layers, each segmented by a 
Icni X 1cm grid. The Hadronic calorimeter has an inner radius of 146cto and has 
thirtyScm thick layers segmented by a 2cm x 2cm grid. The muon calorimeter 
is made of steel-235, has an inner radius of 300 cm and has 22 lOcm thick layers 
segmented by a 10cm, x 10c?n grid. The detector has a 5 Tesla magnetic field. All 
calorimeters are fully sensitive. The cones are made of Tungsten surrounding 
the beamline and form an angle of 15 degrees. See Appendix |A.3| for Higgs 
events simulated using this detector model. 



simulated Event Counts With Total Luminosity = lOOQpb 
Channel: Total Events 
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Figure 2: Simulated event counts for a scan across a 126.0 GeV Higgs peak with 
a 4.2 MeV wide Gaussian beam spread, counting all events except for Z° — > lygiJg 
decays. Data is taken in a 60 MeV range centered on the Higgs mass in bins 
separated by the beam width of 4.2 MeV. Total luminosity is 1 fb~^. Event 
counts are calculated as Poisson-distributed random variables and the data is 
fit to a Breit-Wigner convoluted with a Gaussian peak plus linear background. 
Fitted values of the free parameters are in Table [21 



4.1 Low-Mass Z bosons 

Fortunately, this background is reducible. The s-channel resonance production 
of Higgs bosons only happens with a center of mass energy within a few MeV 
of its peak. Z bosons however are produced in several different processes with 
a wide range of masses, as seen in Figure [3] At an s-channel Higgs factory 
muon collider, Z bosons are primarily produced as real, on-shell bosons along 
with an intial state photon that makes up the difference in energy between the 



Higgs s-channel and the Z mass (Fig. 4b ) . There is also a small number of very 
low mass Z bosons produced in a Drell-Yan process. The only events that are 
theoretically indistinguishable from Higgs events are those where a virtual Z is 
produced at the center of mass energy and decays into a channel shared with 



the Higgs (Fig. 4a I 



Before looking into how the kinematics of these events might differ from 
Higgs events, the simple thing to do is a cut on the total energy potentially 
visible to the detector. This is accomplished by summing the energies of all 
final state particles which pass a cos6' < 0.94 cut and finding the energy cut 
which maximizes S/\/B. The cos 9 cut is effective because most of the high- 
energy initial state radiation is colinear with the beam. We use a cut of Etotai > 
98.0 GeV, which selects 79.2% of the Higgs signal events and 41.9% of the Z 
background. This results in an effective Higgs cross section of 22.4 pb and a 
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Figure 3: Z boson masses in 10,000 PYTHIA-simulated /i^/J.^ — > Z events at 
%fs = VlhXsGeV . The low-mass region is dominated by the Drell-Yan process. 
There is a peak around the Z mass where intial-state Bremsstrahlung radiation 
allows the creation of an on-shcU Z. The third region of interest is the peak at 
125GeV, the center of mass energy. This represents a process with no initial 
state radiation where the off-shell Z's produced are indistinguishable from the 
Higgs. 



background of 126.4 pb, bringing S/\~B to 1.99. Figure Isl shows simulated data 
using these results, with a fitted width of 5.57 ± 1.33 MeV and an error in the 
mass measurement of —0.02 ± 0.14 MeV. This simple cut has already proven 
to be a marginal improvement but there is much more that can be done by 
focusing on individual decay channels. 

4.2 hi 

Table [l] compares the branching ratios and cross sections of the Z background 
with the Higgs signal. The largest Higgs decay channel is H^ — > 66, which makes 
up 58% of Higgs decays at this mass, a branching fraction proportionally large 
to Br{Z^ — > hb) — 15.2%. We assume a b-tagging efficiency and purity of 1, so 
the cross sections for the decays are 16.5 and 57.2 pb, respectively. The fitted 
values for the mass, width and branching ratio of the Higgs using b-tagging are 
shown in Table [2] and a fit to simulated data can be found in Appendix |A.1[ 

In both signal and background the visible energy spectrum is very similar 
to the spectrum of the combined channels, so the same total energy cut of 
Etot > 98.0Gey maximizes S/y/B. Cuts on the event shape, the magnitude of 
the thrust and major axis, can further enhance the signal. The event shape is 
calculated by finding the axis which maximizes the sum of all particle momenta 
projected onto a single axis, called the 'thrust axis'. This is then repeated for 
an axis perpendicular to the first and then a third orthogonal to both. The 
thrust is the normalized sum of the projection of all particle momenta against 
the thrust axis and the major axis value is the normalized sum of the projections 





(a) Irreducible background: 

/i+^~ -5> Zj^* with Mz* = '/a. 



(b) Reducible background: 
/i+/i" -> ^",7 with Mzo < Mfjo. 



Figure 4: Standard Model backgrounds at a /i^/x collider operating at -^s 
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Table 1: Branching fractions and effective cross sections for Standard Model 
decay modes of Higgs and Z bosons. Higgs cross sections are calculated as the 
peak value of the Higgs peak Breit-Wigner convoluted with a Gaussian of width 
4.2 MeV to simulate the effect of beam smearing. Branching fractions are taken 
from PYTHIA 6.4 event generation output. 



against the secondary axis. Because the Higgs is never created in events with 
significant beamstrahlung it is always produced with low momentum. Z bosons 
produced with mass lower than the beam center-of-mass energy are 'boosted' by 
the beamstrahlung photon. This boost lowers the thrust and raises the major 
axis values, so it is a useful indicator for channels with particular event shape 
profiles. 

Figure |6] shows the signal and background thrust and major axes before and 
after cutting on the total energy and event shape values. The cuts were made 
by selecting events with Etot > 98.0GeV, thrust values between 0.94 and 1.0 
and major axis values between 0.0 and 0.20. We continue to assume perfect 
b-tagging. These cuts reduce the bb signal by 52% and the background by 15%, 
bringing the effective cross sections to 8.64 and 8.45 pb respectively. This brings 
the S/^/~B ratio to 2.97, a dramatic improvement over simple energy cuts or b- 
tagging alone. Figure [7| shows a simulated scan of the Higgs peak with a fit to 
a Breit-Wigner convoluted with a Gaussian. 



simulated Event Counts With Total Luminosity = lOOOpb ^ 
Channel: Total Energy Cut 
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Figure 5: Simulated event counts for a scan across a 126.0 GeV Higgs peak 
with a 4.2 MeV wide Gaussian beam spread, counting all events with a total 
energy of at least 98.0 GeV visible to the detector. Data is taken in a 60 MeV 
range centered on the Higgs mass in bins separated by the beam width of 4.2 
MeV. Event counts are calculated as Poisson-distributed random variables and 
the data is fit to a Gaussian peak plus linear background. The fit width is 
5.16 ± 0.24 MeV and the error in the mass measurement is 0.26 ± 0.19 MeV. 



4.3 H^ -^ WW* 

There are several channels with very little physics background that are of im- 
portance, despite their smaller cross sections. One of these is the H^ —>■ WW* 
decay mode, with a branching fraction of 0.226 (cross section 6.39 pb) and no 
real background from the corresponding Z decays. The W boson decays into a 
charged lepton and corresponding neutrino 32.4% of the time, with effectively 
equal rates for each type of lepton. The majority of the remaining branching 
fraction is the decay into pairs of light quarks. While it is certainly possible 
to reconstruct W bosons from four-jet events, in this report we focus on the 
decays with missing energy in the form of neutrinos since they can be identified 
by the presence of one or two isolated leptons and missing energy and are the 
most common. Further study will be required for a detailed analysis of the 
four-jet case. Since the W boson decays into a lepton and neutrino 32.4% of 
the time and we require at least one such decay between a pair of W's, these 
make up 54.3% of WW* events. Thus the theoretical cross section is 6.39 pb 
with virtually no background. 

Because the detector will have a non-sensitive cone, there will be a small 
amount of 'fake' background, eg. when the photon in the decay ^Ji~^ n~ -^ Z'^ + 
J ^ £~^ + i" boosts the two leptons and disappears into the cone as missing 
energy. Figure [26] in Appendix |A.3| shows an example event display for a WW* 
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Figure 6: Effects of event sliape and energy cuts on Higgs bb signal and back- 
ground. Cuts were made by selecting events with total energy Etot > 98.0Gey 
visible to the detector, thrust between 0.94 and 1.0 and major axis between 0.0 
and 0.2. The signal is reduced to 52% and the background to 15%. 



decay into two leptons and illustrates the characteristic missing energy of these 
events. It is difficult to estimate the true background from processes such as 
these, but given the low branching ratios of Z'^ to lepton pairs and the kinematic 
and geometric constraints for 'fake' background, it is safe to assume that the 
background will be fairly low in this channel. Therefore we use the rate assumed 
by Han et al"^], a cross-section of 0.051 pb. Plots of simulated data for the WW* 
channel can be found in Appendix |A.1| and fitted values in Table |3] 

4.4 T~^T~ 

The T'^T^ channel is dominated by the background, but the Higgs branching 
ratio of 0.071 is not insignificant. The Z° — > t'^t~ process has a branching 
ratio of 0.034, giving it an effective cross section of 12.8 pb, compared to the 
2.01 pb cross section for the Higgs. However, the boost given to the lower mass 
Z bosons means the background can be further distinguished using total energy 
and event shape parameters. 

The T is a short-lived particle and every r decay channel involves the pro- 
duction of a T neutrino. This makes the total visible energy less useful as a cut 
parameter than it was for 66, since there are random amounts of missing energy. 
We require at least 60.0 GeV to be visible because background dominates below 
this value due to boosted Z's. Event shape parameters, however, are very useful 
here since t decays typically do not create a widespread shower. We require the 
thrust to be between 0.999 and 1.0 and the major axis to be between 0.007 and 
0.03. This cut reduces the signal to 78% of its original value and the background 
to 39%, bringing the Higgs cross section to 1.58 pb and the background to 4.97 
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Simulated Event Counts With Total Luminosity = lOOOpb 
Channel: bb With Energy and Event Shape Cuts 
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Figure 7: Simulated event counts for a scan across a 126.0 GeV Higgs peak 
with a 4.2 MeV wide Gaussian beam spread, counting X ^ bb events with a 
total energy of at least 98.0 GeV visible to the detector and cutting on event 
shape parameters. Data is taken in a 60 MeV range centered on the Higgs mass 
in bins separated by the beam width of 4.2 MeV. Event counts are calculated 
as Poisson-distributed random variables and the data is fit to a Breit-Wigner 
convoluted with a Gaussian plus linear background. The fit width is 4.78 ±0.48 
MeV, the error in the mass measurement is 0.01 ± 0.05 MeV and the branching 
ratio is measured at 0.271 ± 0.001. Total luminosity is 100Qpb~^, or 71Apb~^ 
per point. 
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Figure 8: Effects of event shape and energy cuts on Higgs t^t^ signal and 
background. Cuts were made by selecting events with total energy Efot > 
GO.OGeV visible to the detector, thrust between 0.999 and 1.0 and major axis 
between 0.07 and 0.032. The signal is reduced to 78% and the background to 
39%. 
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Table 2: Branching fractions, cross sections before and after cuts and S/\/B for 
the channels studied. 



pb, as seen in Figure [8J The cut is specific enough that it is not necessary to 
assume anything else about the events, such as a perfect t+t~ tag. Fewer than 
0.2% of the Higgs decays that pass the cut are not t+t~ events and only 6.4% 
of the background events that pass are misidentified. The effective background 
cross section above is calculated from all the events which pass the cut. Plots 
of simulated data can be found in Appendix | A . 1 1 and fitted values can be found 
in Table H 

4.5 H^ -^ 77 

The final channel examined in this report is the H'^ — ^ 77 channel. The Higgs 
branching fraction for this channel is only 0.3%, but the events can't be easily 
identified by selecting events with two photons with equal energy adding up 
to ^/s and high thrust. About 10% of these events are lost when one or both 
photons hit the cone and there is no background so the cross section is 0.077 pb. 
The high purity of this channel is a great advantage, but the small cross section 
makes it impractical for scanning the beam energy to find the Higgs peak as it 
takes a great deal of luminosity to expect more than a few events on the peak. 
This channel will require much luminosity but may prove very useful for precise 
measurements of the Higgs. 



5 Higgs Measurements 

In the previous section wc fit simulated data to extract the properties of the 
Higgs. While it is clear that the bb and WW* channels will be the most useful 
for measuring Higgs properties, particularly with lower luminosities, the results 
of these fits are not reliable estimates of the achievable accuracy and preci- 
sion of a muon collider. The values quoted were individual samples from trials 
that varied significantly in both accuracy and precision and which used the ap- 
proximation that the background cross section, luminosity per point and beam 
resolution arc well-known parameters. In this section we maintain this assump- 
tion and estimate the achievable accuracy and luminosity dependence of these 
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Table 3: Fitted values of Higgs decay width, mass and branching ratio from 
simulated data. Mass values are the difference between the measured mass 
and the true mass of 126,000 MeV. Total integrated luminosity was 1 fb^^, or 
71Apb~^ per data point. 



measurements . 

5.1 Measurements With the bb Channel 

The uncertainties in the measured values do not always reflect the accuracies of 
the measurements or their statistical variance from experiment to experiment. 
To get a better estimation we repeated the experiment of simulating 1 fb^^ 
of data and fitting it forty times. Figure [9] shows the results of this in box- 
and-whisker plots for a range of integrated luminosities. To reiterate, each 
experiment simulates taking data in a 60 MeV range around the Higgs peak with 
14 bins separated by the beam width of 4.2 MeV. The integrated luminosity is 
the sum of luminosity taken in each bin. 

These plots demonstrate that our simplistic simulation and fitting experi- 
ment is on average accurate, but the statistical variance is high. While a more 
thorough analysis may provide more consistent results, we conclude here that at 
a given luminosity, the Higgs parameters can be measured to within the inner- 
quartile range given. At an integrated luminosity of 1 fb~^, we can use the 
bb channel with energy and event shape cuts to accurately measure the mass 
of the Higgs to within 0.3 MeV, the partial width to within 0.9 MeV and the 
branching ratio to within 0.09. 



5.2 Measurements with the WW* Channel 

We performed the same simulated experiments using our estimated cross sec- 

We find 
channel 



tions for the WW* channel and background, as shown in Figure 10 

that with an integrated luminosity of 1 fb~^, we can use the WW 

with a lepton and missing energy to accurately measure the mass of the Higgs 

to within 0.38 MeV, the partial width to within 0.75 MeV and the branching 

ratio to within 0.02. These values can be found in Table |4] 
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Fitted Values of Higgs Parameters vs. Integrated Luminosity 
66 Events With Event Shape and Energy Cuts 
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Figure 9: Box-and-whisker plots of fitted values of the Higgs mass, bb partial 
width and bb branching ratio for 40 experiments at each luminosity. Integrated 
luminosity is the total luminosity taken in 14 bins 4.2 MeV apart in a 60 MeV 
range centered on the Higgs mass. The boxes extend to the upper and lower 
quartiles of the data and the 'whiskers' extend to the most extreme value within 
1.5 times the inner- quart ile range. 
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Fitted Values of Higgs Parameters vs. Integrated Luminosity 
WW* Events With Lepton and Missing Energy 
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Figure 10: Box-and-whisker plots of fitted values of the Higgs mass, WW* 
partial width, and WW* branching ratio for 40 experiments at each luminosity. 
Integrated luminosity is the total luminosity taken in 14 bins 4.2 MeV apart in 
a 60 MeV range centered on the Higgs mass. The boxes extend to the upper 
and lower quartiles of the data and the 'whiskers' extend to the most extreme 
value within 1.5 times the inner-quartile range. 



16 



Channel 


SMh (McV) 


STh (MeV) 


6Br{H'> - 


-^X) 


bb 


0.30 


0.60 


0.09 


WW* 


0.40 


0.75 


0.02 


Combined 


0.25 


0.45 


— 



Table 4: Accuracy of fitting parameters for simulated Higgs data. Values rep- 
resent the inner quartile range (25% to 75%) of the values of 40 simulated 
experiments using 1 fb^^ total integrated luminosity. The combined values 
were calculated after each experiment using a weighted average. 



5.3 Combining Channels 

To measure the Higgs mass and total width more precisely, we took advantage 
of both channels. We did this by simulating data for both channels at the same 
time and taking their average, weighted by the uncertainty in the fits. For 
example, the formula used for the width was: 



STh = 



ST^bb 



ST,-, + 6r 



-r 



sn 



WW 



WW 



ST., + ST 



WW 



~^bb 



(4) 



As shown in Figure 11 the mass measurement was found to be accurate within 
0.25 MeV and the total width was accurate within 0.45 MeV. All the estimated 
accuracies can be found in Table |4l 



6 Finding the Higgs Peak 

Before precise measurements of the Higgs can be made it will be necessary to 
search for the Higgs by scanning the beam energy and looking for a significant 
signal. Because s-channel resonant production only happens with the beam 
center-of-mass energy directly on the Higgs mass, this is not a trivial matter 
and it is important to know how much luminosity will be required to find the 
Higgs, what channels will be most useful, and what beam momentum spread 
will be ideal; too wide and the effective Higgs cross section will be miniscule, 
too narrow and it will require too many data points. 

To answer this question we devised a simple scanning scheme and calculated 
how much luminosity would be needed to guarantee finding a significant Higgs 
signal. We begin by assuming the Higgs mass has been measured at 126.0 GeV 
to within 100 MeV, based on estimates of the potential resolution possible at 
the LHC.[3] We then calculate how much luminosity will be needed at each 
point to guarantee, to a given confidence level, a statisitcally significant signal 
assuming the beam center of mass energy is exactly at the Higgs mass. Then 
we conduct a search guided by the existing mass measurement, taking data at 
energies calculated to add the most probability to finding the Higgs. 

6.1 Guaranteeing a Significant Signal 

The first step is finding out how much luminosity will be needed to guarantee 
a significant signal if we are taking data on the peak. We will assume we are 
requiring a signal with 5cr significance. To do this, we first calculate N, the 
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Fitted Values of Higgs Parameters vs. Integrated Luminosity 
bb and 'nw* Weighted Averages 
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Figure 11: Box- and- whisker plots of fitted values of the Higgs mass and total 
width for 40 experiments at each luminosity. Integrated luminosity is the total 
luminosity taken in 14 bins 4.2 MeV apart in a 60 MeV range centered on the 
Higgs mass. The boxes extend to the upper and lower quartiles of the data 
and the 'whiskers' extend to the most extreme value within 1.5 times the inner- 
quartile range. 
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number of events that must be observed to count as a significant signal, given a 
certain luminosity and background. We use the Poissonian distribution with an 
expected value of £ x tTf, events to calculate the probability p of a given signal 
and require an N such that p is less than 5cr where ab is the background cross 
section. 



P 

In 



/•OO 

/ Pr(X = n\C X o-fc) dn (5) 

Jn 



Then we calculate the confidence level a, the probability of not seeing more than 
N events given we are taking data on the peak, where the expected number of 
events is £ x (at + (Js) and as is the signal cross section. 



a= Pr(X = n\£ x {at + <Ts)) dn (6) 

"'0 

With the Poissonian, these integrals reduce to a sum. We calculate the lumi- 
nosity required for the confidence level a to reach a range of values (Icr, 1.64(7, 
2a, 2.58cr, 3cr, 4cr, 5a). 

6.2 Search Strategy 

We can now evaluate the probability of finding a significant Higgs signal at a 
given energy and luminosity based on the confidence level a provided by that 
luminosity and a Gaussian probability distribution based on the LHC measure- 
ments. We then order the search by maximizing the probability of finding a 
significant signal at each step, taking enough data at each point to raise a to 
its next value. Data points are separated in energy by the width of the beam 
spread Gaussian to minimize the risk of missing the Higgs peak between data 
points while also minimizing the number of data points needed, since sums of 
Gaussians separated by their standard deviation form an approximately fiat 
distribution. We then calculated how much total integrated luminosity would 
be required for this search strategy to guarantee finding the Higgs peak at a 
range of luminosities, given a certain signal and background and using a range 
of beam widths. 

6.3 Required Luminosity 

First, we looked at the results of using simple event counting, using all events 



except for /i'^/x — >■ .Z'' — )■ vivi. Figure 12 makes it clear that total event count 



ing would not be practical; it would require on the order of inverse femtobarns of 
integrated luminosity to guarantee finding the Higgs at a 3a confidence level. It 
does however reveal that a beam width near the actual Higgs width would offer 
the ideal balance. We then looked at the improvement that could be achieved 
with the energy cuts described above (See Figure [13| . Required luminosities 
for these and the following channels are listed in Table [5J Additional plots of 
luinosity can be found in Appendix ] A. 2 [ 

We next tried the search strategy using tagging of bb events and the effect 
of implementing the energy and event shape cuts described above. As shown 
in Figures [21] and [22j the cuts reduced the required luminosity by a factor of 
about two. As expected, a beam width near the actual Higgs width continues 
to be the most effective. 
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Probability of NOT Finding 5a Higgs Signal by Luminosity: 
Simple Event Counting 




Luminosity - pb 



Figure 12: Probability of not finding a significant Higgs signal using our search 
strategy as a function of luminosity for a range of beam energy resolutions. 
Calculated for a search using simple event counting only. 
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Probability of NOT Finding Scr Higgs Signal by Luminosity: 
Total Events With Energy Cut 
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Figure 13: Probability of not finding a significant Higgs signal using our search 
strategy as a function of luminosity for a range of beam energy resolutions. 
Calculated for a search using simple event counting and requiring events with 
total visible energy Efot > 98GeV. 
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Probability of NOT Finding 5a Higgs Signal by Luminosity 
66 Tagging With Cuts, H^' -^WW* With Lepton and Missing Energy 
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Figure 14: Probability of not finding a significant Higgs signal using our search 
strategy as a function of luminosity for a range of beam energy resolutions. 
Calculated for a search combining the bb and WW* channels, using the cuts 
described above. 



Finally, we looked at H^ — >■ WW* events with a lepton and missing energy, 
as described above. The high purity and relatively large cross section of this 
channel makes it significantly better than other channels for searching for the 
Higgs (Figure 23). Interestingly, it also appears that the low background re- 
moves the necessity for a beam with a width near the Higgs decay width; beam 
widths in the range from 4.2 MeV to 25 MeV fare similarly with this search 
strategy in terms of luminosity. 



6.4 Combined Channels 

Fortunately, there is no need to use only one channel or another when searching 
for the Higgs, so next we examine the effects of using both the bb and WW* 
channels in our search strategy. As seen in Figure [T4| this offers improvement 
over using any one channel alone. This is because the probability of the back- 
ground fluctuating significantly in two channels at once is equal to the product 
of the probabilities for each channel fluctuating on their own. Again we see that 
the beam resolution makes little difference and while wider beam spreads seem 
to require slightly less luminosity, narrower beams will measure the Higgs peak 
with greater precision, reducing luminosity needed for more precise measure- 
ments. 
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Channel 


CTsig (pb) 
CTbkgr (pb) 


Luminosity Required {pb ^) 


Ci = 3(T 


CL = 5a 


Total 


as ^ 28.3 
ah = 301.4 


1,723 


3,840 


Total (Cut) 


as = 22.4 
ah = 126.4 


1,193 


2,666 


bb 


(Ts — 16.5 
ah = 57.2 


1,033 


2,317 


bb (Cut) 


(Ts = 8.64 

ah = 8.45 


697 


1,593 


WW* 


(Ts = 6.39 
(T,, = 0.05 


146 


389 


WW* (Cut) 


(Ts = 3.35 
(Tt = 0.05 


325 


812 


bb, WW* 


CTs = 

CTb = 


105 


368 



Table 5: Required luminosity to guarantee finding a 5a Higgs signal with con- 
fidence level a — 3cr, 5a. 

7 Discussion and Conclusion 

The Higgs boson is a particle of fundamental importance to physics and measur- 
ing its properties with precision will allow us to probe the limits of the Standard 
Model and may point the way towards non-Standard model physics. Using sim- 
ple estimates of physics backgrounds and separable signal we have estimated 
that with 1 fb~^ of integrated luminosity a hypothetical muon collider Higgs 
factory operating at the Higgs s-channel resonance could measure the mass of a 
Standard Model 126 GeV Higgs to within 0.25 MeV and its total width to within 
0.45 MeV. We estimated that with a beam spread of 4.2 MeV, approximately 
368 pb^^ total integrated luminosity would be required to guarantee locating the 
narrow Higgs peak. We believe that these preliminary results strongly motivate 
further research and development towards the construction of a muon collider 
Higgs factory. 

Our estimations assume that there is no machine-induced background and 
that the detector has excellent tracking and calorimetry. Our results demon- 
strate the value of the high Higgs cross section and narrow beam energy spread 
available at a muon collider. These two factors enable the direct measurement 
of the Higgs mass and width by scanning the Higgs s-channel resonance, which 
is not possible at any e+e^ collider. Our study of the physics-induced back- 
ground and separation of the Higgs signal showed that significant reduction of 
the physics background can be achieved by a detector with high energy and 
spatial resolution. We believe that this report justifies more in-depth analy- 
sis of Higgs channels and their backgrounds, for example the reconstruction 
of iJ" -H- WW* — >■ 4j events using learning algorithms or the application of 
flavor-tagging techniques to tag bb events. 

Machine-induced backgrounds, mainly from muon decays in the beam, present 
an additional difficulty which has not yet been studied in great detail. We be- 
lieve that in addition to significant shielding in the detector cone and endcaps. 
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it may be important to have a calorimeter with high spatial and temporal reso- 
lution. Our results motivate an in-depth analysis of the machine-induced back- 
ground including simulation in a highly segmented, totally-active, dual readout 
calorimeter such as the MCDRCalOl detector concept. 
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A Additional Figures 
A.l Simulated Event Counts 



simulated Event Counts With Total Luminosity = lOOOpb 
Ctiannel: 66 Tagging 
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Figure 15: Simulated event counts for a scan across a 126.0 GeV Higgs peak 
with a 4.2 MeV wide Gaussian beam spread, counting all X —i' bb events. Data 
is taken in a 60 MeV range centered on the Higgs mass in bins separated by 
the beam width of 4.2 MeV. Event counts are calculated as Poisson-distributed 
random variables and the data is fit to a Breit-Wigner convoluted with a Gaus- 
sian plus linear background. The fit width is 3.49 ± 1.83 MeV, the error in the 
mass measurement is —0.06 ±0.19 MeV and the branching ratio is measured at 
0.67 ±0.05. Total luminosity is 1000p5"\ or 71Apb~^ per point. 



24 



450 
400 
350 
300 

a 250 

S200 

150 
100 
50 



Simulated Event Counts With Total Luminosity = lOOOpfc 
Channel: HW 
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Figure 16: Simulated event counts for a scan across a 126.0 GeV Higgs peak 
with a 4.2 MeV wide Gaussian beam spread, counting all H^ — >■ WW* events 
with a minimal background. Data is taken in a 60 MeV range centered on the 
Higgs mass in bins separated by the beam width of 4.2 MeV. Event counts 
are calculated as Poisson-distributed random variables and the data is fit to a 
Breit-Wigner convoluted with a Gaussian plus linear background. The fit width 
is 4.06 ± 0.24 MeV, the error in the mass measurement is 0.00 ± 0.07 MeV and 
the branching ratio is measured at 0.217 ± 0.001. Total luminosity is 1000pb~^, 
or 71Apb^^ per point. 



simulated Event Counts With Total Luminosity = lOOOpfc 
Channel: WTr* With Lepton and Missing Energy 
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Figure 17: Simulated event counts for a scan across a 126.0 GeV Higgs peak 
with a 4.2 MeV wide Gaussian beam spread, counting all H'^ — >■ WW* — >■ 
lepton + missing energy events with a minimal background. Data is taken in 
a 60 MeV range centered on the Higgs mass in bins separated by the beam 
width of 4.2 MeV. Event counts are calculated as Poisson-distributed random 
variables and the data is fit to a Breit-Wigner convoluted with a Gaussian 
plus linear background. The fit width is 3.96 ± 0.17 MeV, the error in the 
mass measurement is —0.16 ± 0.04 MeV and the branching ratio is measured at 
0.1271 ± 0.0002. Total luminosity is 1000p6~\ or 71 Apb'^ per point. 
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simulated Event Counts With Total Luminosity = lOOOpb ' 
Channel: r^r" 
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Figure 18: Simulated event counts for a scan across a 126.0 GeV Higgs peak 
with a 4.2 MeV wide Gaussian beam spread, counting all X — )■ t^t~ events. 
Data is taken in a 60 MeV range centered on the Higgs mass in bins separated by 
the beam width of 4.2 MeV. Event counts are calculated as Poisson-distributed 
random variables and the data is fit to a Breit-Wigner convoluted with a Gaus- 
sian plus linear background. The fit width is 4.82 ± 4.46 MeV, the error in the 
mass measurement is — 0.54 ± 0.47 MeV and the branching ratio is measured at 
0.0623 ± 0.0005. Total luminosity is lOOOpb''^, or 71 Apb'^ per point. 



simulated Event Counts With Total Luminosity = lOOQpb ^ 
Channel: t^ t^ With Energy and Event Shape Cuts 
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Figure 19: Simulated event counts for a scan across a 126.0 GeV Higgs peak 
with a 4.2 MeV wide Gaussian beam spread, counting X — >■ t^t~ events with 
a total energy of at least 60.0 GeV visible to the detector and cutting on event 
shape parameters. Data is taken in a 60 MeV range centered on the Higgs mass 
in bins separated by the beam width of 4.2 MeV. Event counts are calculated 
as Poisson-distributed random variables and the data is fit to a Breit-Wigner 
convoluted with a Gaussian plus linear background. The fit width is 0.84 ±2.97 
MeV, the error in the mass measurement is 1.07 ±0.30 MeV and the branching 
ratio is measured at 0.24 ± 0.23. Total luminosity is lOOOpb^^, or 71.4p6~^ per 
point. 
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simulated Event Counts With Total Luminosity = lOOOpb ' 
Channel: -,-) 




Figure 20: Simulated event counts for a scan across a 126.0 GeV Higgs peak 
with a 4.2 MeV wide Gaussian beam spread, counting all i?° — >■ 77 events. Data 
is taken in a 60 MeV range centered on the Higgs mass in bins separated by 
the beam width of 4.2 MeV. Event counts are calculated as Poisson-distributed 
random variables and the data is fit to a Breit-Wigner convoluted with a Gaus- 
sian plus linear background. The fit width is 2.85 ± 5.73 MeV, the error in the 
mass measurement is —0.6 ± 0.9 MeV and the branching ratio is measured at 
0.0035 ± 0.0001. Total luminosity is lOOOpb'^, or 71 Apb-^ per point. 

A. 2 Search Strategy 



Probability of NOT Finding 5a Higgs Signal by Luminosity: 
Perfect bb Event Tagging 
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Figure 21: Probability of not finding a significant Higgs signal using our search 
strategy as a function of luminosity for a range of beam energy resolutions. 
Calculated for a search using perfect bb tagging. 
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Probability of NOT Finding Su Higgs Signal by Luminosity: 
Perfect bb Tagging, Event Shape and Energy Cut 
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Figure 22: Probability of not finding a significant Higgs signal using our search 
strategy as a function of luminosity for a range of beam energy resolutions. 
Calculated for a search using perfect bb tagging and event shape and energy 
cuts. 



Probability of NOT Finding 5(t Higgs Signal by Luminosity: 
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Figure 23: Probability of not finding a significant Higgs signal using our search 
strategy as a function of luminosity for a range of beam energy resolutions. 
Calculated for a search using i/° -^ WW* -> lepton + missing energy events 
with minimal background. 
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A. 3 Event Displays 

The following event displays were generated using the WIRED4 plug-in for JAS3 
using fi+fj,- -^ H° events with Vs = Mh = 125.0GeV generated by PYTHIA 
6.4 and simulated in the MCDRCalOl concept detector using SLIC. 




Figure 24: Example _ff° — > bb event display. 
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Figure 25: Example if'^ — >■ bb event display with fisheye projection to illustrate 
the displaced b vertices. 




Figure 26: Example H° — > WW* -> e + i^e + t+ + t/^. event display. Yellow 
lines indicate neutrinos. 
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Figure 27: Example i/° -^ t^t event display. 
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